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ABSTRACT 


A simple theoretical model explaining basic properties of radar Imaging of 
underwiter bottom topography in tidal channels i9 presented. The surface 
mgliness modulation Is described by weak hydrodynamic interaction theory in the 
relaxation time approximation. In contrast to previous theories on short wave 
modulation by long ocean waves, a different approximation has to be used when 
describing short wave modulation by tidal flow over underwater bottom topography. 
The modulation depth is in this case proportional to the relaxation time of the 
Brigg waves. The large modulation of radar reflectivity observed in SEASAT-SAR 
imagery of sand banks in the Southern Bight of the North Sea can be explained by 
issumLng that the relaxation time of 34 cm Bragg waves is of the order of 30-40 
seconds . 


1. INTRODUCTION 


It has been known for more than 10 years that real aperture radar (RAR> imagery 
taken over sea areas with strong tidal currents (tidal channels) sometimes shows 
features that seem to be related to underwater bottom topography (de Loor and van 
Huiten, 1978; de Loor, 1981; Mcleish et. al., 1981). The same phenomenon has also 
been observed in synthetic aperture radar (SAR) imagery obtained by the SEASAT 
satellite (Fu and Holt, 1982; Lodge, 1983; Kenyon, 1983; Lyzenga et. al., 1983). 
The sea-floor topography (bathymetry) causing these radar signatures sometimes 
lies tens of meters below the sea surface. An example of such Imagery Is the 
SEaSAT-SAR image shown In Fig. 1. 
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PIG.l: Digitally processed SEASAT-SAR laage of the Southern Bight of the North 
Sea froa orbit 762 (Aug. 19, 1978, 6:46 UT) with fraae center at 51 # X9* N, I^'e. 
The land area In lower left hand corner la the English coast near Raasgate. The 
V-shaped feature In the center are the sand banks South Palls and Sandettle. 
South Palls U about 30 ka long, 600-800 a broad and rises froa the sea floor of 
a depth of about 30-40 a to within 7 a of the sea surface. 
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The imaging of underwater bottom topogr ipny by real or synthetic aperture radar 
seems not to be understandable at first sight, because the penetration depth of 
the elect romagnet lc waves emitted by the radar Into sea water is only of the 
order ot millimeters to centimeters. Therefore, surface effects related to the 
sea-floor topography must be responsible for the radir Imaging. 

Since radar signatures of bathymetry are only observed when strong (tidal) 
currents are present, we are led to the hypothesis thit these radar signatures 
are caused by surface current variations associated with underwater bottom 
topography. Variable surface currents modify the short-scale surface roughness, 
and this roughness variation Is sensed by the radar. Since most Imaging radars 
operate at Incidence angles between 20 and 70 degrees the radar backscatterlng Ls 
dominated by Bragg scattering (Valenzuela, 1978). Consequentely , croes-section 
modulation results from the modulation of the spectral energy of the Bragg waves. 

Cross section modulation determines the Image Intensity or grey tone level 
variations in real iperLure radar Imagery. However, in case of synthetic aperture 
radar. In addition to this amplitude modulation, also phase modulation or 
velocity bunching contributes to the imaging (Larson et.al. 1976; Alpers an! 
Rufenach, 1979; Alpers et.ii. 1981). In most cases of SAR imaging of bathymetry, 
however, the contribution ot velocity bunching to the ImagLng mechanism Is small 
(Alpers and Hennings, 1984). Therefore we will discuss Ln thLs paper only 
amplitude or cross sectLon modulation. 


2. SATHYMSfRY - CURRENT INTERACTION 

The Interaction of a 3-dimenslonai time-variable current field with 3-dimensional 
underwater bottom topography (bathymetry) can sometimes be a very complex 
process which does not allow a simple mathematical description. 

Nevertheless, in this paper we make the simplest possible assumption that the 
current flow over the bathymetry is laminar, free of any vertical current shear 
and only weakly time dependent. Furthermore, we assume that the tidal velocity 

component normal to the direction of the underwater bank obeys the continuity 
equation 
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and that the parallel component remains constant, 


= const 


Here d(x^) denotes a depth profile along a line perpendicular to the ridge 
direction (x^). 

Current measurements or* and off the sand bank South Falls by Venn and Olier 
(1983) in the Southern Bight of the North Sea (north-eastern approach to the 
English Channel), where the tidal current flows across the bank at an oblique 
angle, have confirmed that equations (la) and (lb) are acceptable first order 
approximations. (However, deviations from these simple relations are observed, 
and a more refined theory should account for them). 


3. CURRENT - BRAGG WAVE INTERACTION 

Since the variation of the surface current due to interaction with bathymetry has 
space and time scales that, in general, are small compared with the space and 
time scales of the Bragg waves, the current -Bragg wave Interaction can be 
described by a WKB (Wentzel - Kramers - Brillouin)-type interaction theory. In 
this theory the transport equation, which describes the variation of the spectral 
energy density of short waves in a slowly varying current field, is the action 
balance or radiation balance equation (Hasselmann et al., 1973; Keller and 
Wright, 1975; Alpers and Hasselmann, 1978; Wright, 1978)# This equation reads 


l|t + ^ -II- N_ s ( -- fc) 


where 


- I N (x£/t) = E (x,k,t) / cj' < 3 > 

Is the actton spectrum, E(x,k,t) the wave spectrum, Aj'the intrinsic frequency of 
the wave In a reference system which is locally at rest, x the space variable, k 
the wavenumber and S(x,k,t) a source function. The waves propagate along 
| trajectories in 4-dimensional phase space which are given by the ray equations 

* /.A£ 


ft) 


..r* « - 


'’■***< 


* ^ 





? 


X = 


where 


9cO 

TT 


k - - 


" 9*0 

91 


fO(x,k,t) - 0 ; (k) + k • U (x. 


t) 


ORIGINAL PAGE f? 
OF POOR QUALITY 


(4 ) 
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denotes the wave frequency in the moving medium with variable velocity U(x,t), 


J ; 
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We assume that the variable surface current only leads to small deviations of the 
action density from equilibrium. Therefore, we write the action density N and the 
surface current U as sums of a constant equilibrium term and a time dependent 
perturbation term 


N (x,k,t) = N (k) + ^ N (x,k,t) 


( 6 ) 
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Furthermore, we approximate the source term S by a diagonal operator 


s (x,k,t) = ja £■ 


O N (x,k,t) 
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is called relaxation rate 


where A Is a parameter with dimension (time) 

, -l 

and TTp* M the relaxation time. Physically, L^is a system constant describing 
the response of the wave system to small deviations from equilibrium caused by 
surface current variations. It is determined by the combined effect of wind 
excitation, energy transfer to other waves due to conservative resonant wave - 
wave interaction, and energy loss due to dissipative processes like wave 
breaking. No measurements of the relaxation time In the open ocean exist. 
However, from theory we expect that Is of the order of 10 to 100 wave periods. 
Applied to SEASAT-SAR Brag3 waves, which have a wavelength of 34 cm and a wave 
period of 0.47 s, this means that T* should lie In the range between 4.7 and 47 

seconds. In this paper we consider (or JU. ) as a free parameter. Inserting eqns. 
(6), (7) and (8) into eqn. (2) and only keeping first order terms yields 
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The ne scales or the three terms on the left hand side are given by the local 
timeT* , the acvection time T^and the relaxation time . The local time is of 
the order of the period of the semi-diurnal tide divided by 2 7T, which is 
i 2 . 5/ 2 IT h£?2h, and the advection time i3 given by 


^ = I <Sq 
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where K Is the wavenumber of the bottom topography. 


For tidal flow over large-scale underwater bottom topography, such as sand banks, 
both, the local time T and the advection time , are usually much large than the 
relaxation time . Therefore the first two terms on the left hand side of eqn. 
(9) can be neglected in comparison with the third term. 


This approximation is applicable for discribing short surface wave modulation by 
tidal flow over bathymetry, but it is completely different from the approximation 
applicable for discribing short wave modulation by long surface waves. In this 
case, the local time is given by jr , where SI is the radian frequency of the 
long surface wave. 12 is typically of the order of 1 or 2 seconds, which is small 
compared to the relaxation time . Therefore, the first term in eqn. (9) is 
dominant in this case. 


We now insert eqn. (/f) and the dispersion relation for short surface waves into 
eqn. (9) and consider only those waves which travel towards or away from the 
antenna look direction (Bragg waves). If we define the projection of the antenna 
axis onto the horizontal plane as x-direction, then in case of tidal flow over 
large-scale bottom topography, we obtain the result 


£e = _ 4 "3 u x ( - } 

E 0 f ^ Y 

and in case of long surface waves 

? E _ _ 4 +£ 3 U X 

E o S2. Q x 


(ii) 


( 12 ) 


Hereof denotes the ratio between the group and phase velocity of the short waves. 
For gravity waves we obtain ^ ■ 0.5^ and for capillary waves ^ “ 1.5. In deriving 
eqns. (11) and (12) we have assumed that the spectral energy of the short waves 
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Note, that according to Bragg scattering theory, 

W J_E_ 

*0 E . ' (u) 
where*"'* denotes the relative cross section modulation. 

Comparison of eqns. (11) and (12) shows that in the first case the modulation is 

proportional to and in the second case proportional to » (2 IT) , 

where is the period of the long surface wave. We will show below that for the 
SEASAT-SAR Bragg waves is of the order of 30 to 40 seconds. This means that for 
the same current gradient, the modulation by tidal flow over bathymetry is a 
factor 15-40 larger. On the other hand, the surface current gradient generated 

by tidal flow over bathymetry is usually much smaller than the surface current 

gradient generated by the orbital at Ion of long surface wavt (typically one 
order of magnitude smaller). Therefore we obtain the net result, that the cross 
section modulation ("modulation depth") associated with tidal flow over sand 
banks is of the same order of magnitude, or even larger than the modulation 
associated with long surface waves. 

Physically, this result is a consequence of the fact that In the case of tidal 
flow over bathymetry the weak interaction can act sufficiently long on the short 
wave system. The limit is given by the relaxation time. In case of short wave 
modulation by long surface waves, the short waves have not sufficient time to 
build up a strong modulation because their dweil-time In flow regions with 
positive or negative velocity gradients is determined by the period of the long 
surface waves. 

If the partial time derivative and advectton terms in eqn. (9) are neglected, 
then eqn. (9) states, that the relaxation of the action density perturbation is 
balanced by "straining" exerted on the wave system by the spatially varying 
current JJ. The right hand side of eqn. (9) represents a refraction term. It 
originates from the fact that the spatially variable current JU refracts the short 
waves, i.e., changes their wavenumber. A wavenumber change causes a local 
perturbation of the equilibrium action (or energy) density spectrum of the short 


I 
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waves, because the short wave spectrum varies as a functLon of wavenumber* Thus, 
the steeper N is, as a function of wavenumber, the larger is the modulation. 
Note also that the modulation is proportional to the current gradient. 

In those cases where the advection time is not small compared to the 
relaxation time T , the second term on the left hand side of eqn. (9) has to be 
retained. The inclusion of this term results in the id' 4 * tier? of a low-pass filter 
to the imaging process. This is discussed in detail in Alpers and Hennings (1984), 


4. RELATIONSHIP BETWEEN BATHYMETRY AND CROSS SECTION MODULATION 


By combining the results of sections 2 and 3 we can derive an expression relating 
the cross section modulation to the large scale bottom topography and the tidal 
current field. Inserting eqns. (la), (lb) and (11) into eqn. (11) yields 


KK 008 ^ 005 ^ 


(14) 


Here denotes the angle between the flight and the sand bank direction, the 
angle between the x (| -and the (undisturbed) flow direction, graded the gradient of 
the depth profile in direction perpendicular to the bank crest, d^ the water 
depth outside the bank area, and ju^J the modulus of the undisturbed current 
velocity. 


Inspection of eqn. (14) shows that radar signatures of sand banks always have 
double sign, which means that the radar image is composed of image elements 
having both enhanced and reduced grey levels relative to the local mean. The sign 
of the modulation is such that increased radar reflectivity occurs always on the 
downstream side, and reduced radar reflecting on the upstream side of the sand 
bank. 


Furthermore, eqn. (14) shows that the modulation depth increases with tidal 
velocity and decreases with water depth. The modulation pattern is not correlated 
with the depth profile d, but with d graded. 

In particular, eqn. (14) predicts that the cross section modulation vanishes when 
- the tidal velocity is zero ^|U 0 ( » 0) 
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the bank crest is aligned parallel to the current direction ( = 90°) 

- the band crest is aligned parallel to the radar look direction (4> 
(However, we do not expect that the last statement fully holds, because it is a 
consequence of the fact that our hydrodynamic interaction model does not include 
the interaction between short surface waves travelling in different directions^ 

According to eqn. (14) the modulation depth depends on the relaxation 
time 6 = A * We expect that, to first order, ZT is independent of wind direction. 
However, it is likely that decreases with wind speed. 

In order to obtain an estimate of the relaxation time we analyzed digitally 
processed SEASAT-SAR imagery and made several image intensity scans across sand 
banks in the Southern Bight of the North Sea (north eastern approach to the 
English Channel) and compared them with bathymetry. One scan was made across 
South Falls along the line shown in Fig.l. The arrow indicates the direction of 

the tidal flow at the time of the overflight. The tidal current velocity was 

o 

0.60 c/s and the wind was blowing from 135 N at 4 m/s. This image was processed 
on the MDA (McDonald Dettweiler and Associates) digital 3EASAT processor by the 
Deutsche Forschungs- und Versuchsans talt fur Luft- und Raumfahrt at 
Oberpfaf fenhof en, FRG. The resolution is 25m * 25m (4 looks). 

Our analysis shows that we can explain the measured modulation depth if we assume 
a relaxation time of the order of 30-40 seconds for the SEASAT-SAR Bragg 
waves (for a wind speed of 4 m/s). This value correspoi ds to 60-80 wave periods 
of 34 cm waves. 


$. CONCLUSIONS 

We have confronted the predictions of our imaging theory with existing 
experimental data. It seems that the predictions are largely confirmed. However, 
it should be stressed that we consider this theory only to be a first-order 
theory* 

Finally, we want to suggest the use of radar imagery of bathymetry In tidal 
channels as a practical means for measuring the relaxation time of short surface 
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waves in the open ocean. 
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